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OPTIMAL AND SUBOPTIMAL CONTROL TECHNIQUE FOR 
AIRCRAFT SPIN RECOVERY 

By John W. Young 
Langley Research Center 

SUMMARY 

An analytic investigation has been made of procedures for effecting recovery from 
equilibrium spin conditions for three assumed aircraft configurations. Three approaches 
which utilize conventional aerodynamic controls are investigated. Included are a con- 
stant control recovery mode, optimal recoveries, and a suboptimal control logic pat- 
terned after optimal recovery results. The optimal and suboptimal techniques are shown 
to yield a significant improvement in recovery performance over that attained by using a 
constant control recovery procedure. 

INTRODUCTION 

Many modern fighter aircraft have inertia and aerodynamic characteristics which 
produce unsatisfactory handling qualities during maneuvers at angles of attack near the 
stall region (refs. 1 and 2). Thus, inadvertent entries into the poststall- spin region 
occur, and, unless correct recovery techniques are promptly initiated, steady developed 
spins may take place. Should a fast, flat spin occur, recovery using conventional 
methods is difficult, if not impossible. 

The objective of the present study is to evaluate analytically the effectiveness of 
various procedures for achieving recovery from equilibrium spin conditions. Three 
separate techniques which utilize conventional aerodynamic controls are investigated. 
Included are recoveries for which constant control surface deflections are used, recov- 
eries obtained by using a six-degree-of-freedom optimization program, and recoveries 
using a suboptimal control logic patterned after optimal recovery results. Comparisons 
between the performance of each recovery procedure are given for three different 
assumed aircraft configurations. Recovery performance is evaluated with respect to 
turns, altitude loss, and time required to achieve recovery. 


SYMBOLS 


All aerodynamic data and flight motions are referenced to the body system of axes 
shown in figure 1. 

b wing span 

Cj rolling-moment coefficient 

C m pitching -moment coefficient 

C n yawing-moment coefficient 

C x longitudinal-force coefficient 

Cy side-force coefficient 

C z vertical-force coefficient 

c mean aerodynamic chord 

g acceleration due to gravity 

h altitude 

Ah altitude loss during recovery 

*X’ *Y’ *Z’*XZ body-axis moments and product of inertia about the center of mass 

m mass of airplane 

p,q,r body-axis rolling, pitching, and yawing angular rates (fig. 1) 
q dynamic pressure, ipV 2 

£t 

r L yaw rate limit; constant used in stability augmentation logic (table III) 
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5 wing area 

u, v, w components of airplane resultant velocity along X^, Y b , and Z b body 

axes (fig. 1) 

V resultant velocity of airplane 

X b ,Y b’ Z b body-axis reference system 

a angle of attack 

constant used in stability augmentation logic (table III) 

ojp trim angle of attack; used as a constant in stability augmentation logic 

(table III) 

/3 angle of sideslip 

6 aileron deflection, positive when right aileron trailing edge is down (fig. 1) 

ct 

6 elevator deflection, positive when trailing edge is down (fig. 1) 

V 

6 r rudder deflection, positive when trailing edge is left (fig. 1) 

Euler angles 


p air density 

oj angular rate about center of mass, yP^ + $ + 


Stability derivatives: 




8 6 e 



8 C n 
8 6 r 
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c, -i£l 

*6a 96a 


c x 5e = 


8Cx 

9 5 e 



C Yp 


ac Y 




9 £i 

96 r 


Cm 


6 e 


9 Cm 
3 6 e 


Cm q 


ac m 




ac n 



Cn r - 


ac n 




9 Cn 

96 a 


C Y r 


8Cy 





8 Cy 

"867 


9Cy 

CY «r = 


Cz «e 


££z 

36 e 


A dot over a symbol indicates a derivative with respect to time. 

METHOD 

The effectiveness of various spin recovery techniques was evaluated analytically 
for three different aircraft configurations. The nonlinear aircraft dynamics was repre- 
sented by use of a six-degree-of -freedom model, and the associated equations of motion 
are given in the appendix. 
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Debcription of Aircraft Configurations 

Mass and dimensional characteristics for the assumed aircraft, herein referred 
to as configurations A, B, and C, are shown in table I. Also presented in table I are the 
control surface deflection limits for the respective aircraft. These data were taken 
from references 3 and 4. Configuration A represents a twin-jet swept -wing fighter, and 
configurations B and C represent two different delta-wing fighters. 

Aerodynamic data for configuration A were taken from references 5 and 6. Two 
assumptions were made relating to these data. As shown in reference 5, the rudder and 
pitch-down elevator effectiveness at high angles of attack is so small in magnitude and 
so uncertain in sign that it is practically indistinguishable from zero. Therefore, for 
angles of attack greater than 50°, rudder effectiveness and pitch-down elevator effective- 
ness were assumed to be 0°. Aerodynamic data for configurations B and C were taken 
from references 4 and 7, respectively. 

Spin Characteristics of Aircraft Configurations 

Previous analytical studies (for example, ref. 3) have shown that several spin 
modes may exist for aircraft of the type given in table I. These spin modes can be 
characterized as steep, intermediate, or flat, depending on their angle of attack (the 
higher the angle of attack, the flatter the spin). Spins may be either oscillatory or 
steady in nature. By a steady spin it is meant that the aircraft is in the poststall region 
of its flight envelope and is descending along a helical path, the axis of which is vertical. 
Also, the axis of the helix, its radius, and the descent speed are constant, atmospheric 
density variations with altitude and other external disturbances being neglected. Under 
these conditions the linear and angular velocities are constant in the body set of axes 
(a=/3 = V = p= q = r = 0). 

Since recoveries from flat, steady spins are generally the most difficult to achieve, 
only the flat mode is considered in the analysis. Those spin modes considered are sum- 
marized in table II. As shown therein, two flat spins exist for configuration A, whereas 
only one flat spin exists for configurations B and C. The spin characteristics given are 
either taken from reference 3 or are calculated by using the analytic spin prediction 
method of reference 3. As also shown in table II, only erect spins (positive angle of 
attack) are considered. 


Recovery Techniques 

Three separate recovery techniques which utilize conventional aerodynamic con- 
trols were analyzed. These techniques are referred to as (1) constant control recov- 
eries, (2) optimal recoveries, and (3) suboptimal recoveries. A brief description of 
these recovery techniques is given in this section. 
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The assumed aircraft configurations of table I are all loaded heavily along the 
fuselage (lx - Iy negative). The generally recommended recovery technique (ref. 2) 
for this type of aircraft is to apply aileron with the spin (bank right in a right spin) , 
rudder against the spin, and pitch-up elevator. The constant control recovery mode of 
the current analysis is similar to the generally recommended recovery procedure and 
was included in order to establish a standard for comparison of other recovery methods. 
Recoveries were initiated by using full aileron with the spin and full rudder against the 
spin. Elevator position was determined by computing recoveries at various elevator 
settings and then using that setting which gave the best recovery performance with 
regard to recovery turns, altitude loss, and time. These control deflections were main- 
tained until angle of attack and spin rate were reduced sufficiently to return the vehicle 
to a controllable condition. A stability augmentation or rate-damper control logic was 
then used to damp out oscillations in pitch, roll, and yaw and to trim the aircraft at a 
desired angle of attack below the stall region. The augmentation logic used in the anal- 
ysis is outlined in table III. The significance of various parameters in the table is 
discussed in later sections of the report. 

The second recovery technique involved an application of optimal control theory. 

A computer program of three-degree-of-freedom trajectory optimization (ref. 8) was 
modified to include a full six-degree-of -freedom capability. The computer program was 
used to determine optimal aileron, rudder, and elevator control histories for effecting 
recovery. The program utilized an iterative, gradient approach to determine optimal 
control histories so that some functions of the state variables were minimized while con- 
straints on other state and control variables were met. 

The third recovery technique employed a suboptimal control logic in conjunction 
with the stability augmentation system of table III. During the high angle-of-attack 
phase of a recovery, control histories were patterned after optimal recovery results, 
whereas, during the terminal phase of a recovery, the damper logic of table m was used. 

It should be noticed that the stability augmentation logic of table III was used only 
for constant control and suboptimal recoveries. For optimal recoveries, the optimiza- 
tion technique determined all the control histories. Also, as shown in table III, the 
augmentation system provided damping only about one aircraft axis at a time with no 
control cross connects; that is, the elevator was used only for pitch damping, the ail- 
eron for roll damping, and the rudder for yaw damping. 

RESULTS 

Recovery performance from the equilibrium spin conditions of table II is evaluated 
for the three aircraft configurations of table I. The results obtained by using a constant 
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control, optimal, and suboptimal recovery logic are compared with respect to recovery 
time, turns, and altitude loss. 


Configuration A 

Constant control recovery. - A constant control recovery time history for configu- 
ration A is shown in figure 2. The initial conditions for this recovery are the same as 
those of run 1, table II. As previously outlined, constant control recoveries are made 
by using those combinations of constant control deflections which give the best recovery 
performance. These control settings are maintained until the angle of attack and spin 
rate are reduced sufficiently for a stability augmentation system to damp out pitch, roll, 
and yaw rates and trim the aircraft at a desired angle of attack. (Recovery turns, time, 
and altitude loss refer to the values attained at this point and not that required to regain 
straight and level flight.) 

For the recovery data of figure 2, the best results were obtained by using full rudder 
opposing the spin, full aileron with the spin, and a pitch-up elevator of -10°. As shown 
in figure 2, these control settings were maintained for about 18 sec, at which time the 
stability augmentation or rate-damper logic of table III was activated. 

The augmentation logic is activated for a < o^. The analysis showed that values 
for between 40° and 50° were adequate for all configurations of the study. Although 
this parameter is not overly critical, it must be low enough for the controls to be effec- 
tive. A value of a l = 50° was used for all configurations since this is the angle of 
attack at which the rudder generally becomes effective. 

As is shown in figure 2, for a > a L , 6 r = 0°. This occurs since the rudder was 
assumed to be ineffective for a > 50° and was set at zero for computational purposes. 

In actual practice, the rudder could have been maintained at 6 r = 30° until the yaw rate 
was reduced to zero. 

As is shown in table III, two additional parameters, r L and a T , are required to 
define the stability augmentation system. The yaw-rate limit r L was used to insure 
that the aircraft was out of a spin before an attempt was made to trim the vehicle at a 
desired angle of attack a T . As was the case with « L , recovery performance was not 
found to be overly sensitive to the assumed value for r^, so that values of approximately 
0.2 to 0.6 rad/sec were adequate for all aircraft configurations considered. The value 
chosen for r L was 0.4 rad/sec for all configurations except configuration B. Better 
recovery performance was achieved with configuration B by removing the r L require- 
ment; that is, for a < 50°, the pitch damper was not used and the trim angle of attack 
was commanded. 

In figure 2 the pitch-down elevator deflection at approximately 25 sec is a result 
of the switching of the damper logic to the trim angle -of -attack phase. The desired trim 
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angle of attack otj, was set at 10° in the current study, as this value is well below the 
stall for all aircraft considered. This angle of attack appears to be a reasonable point 
from which to initiate the final pullout following recovery since it resulted in a low-g 
pullup maneuver for all assumed configurations (0.5g to 1.5g). In actual practice, a trim 
angle of attack near the maximum lift condition might be desirable because it would yield 
a more rapid pullup. 

It should be emphasized that the previously described recovery control logic and 
damper logic were not "tuned” to each configuration, since the objective was to investi- 
gate basic principles of recovery control which were not sensitive to configurations. In 
any real control system design, variables such as o^, r L wou ^ probably 

be adjusted for each configuration. In addition, the damper logic for a practical system 
would have to include realistic limits on control surface deflection rates. The near 
"bang -bang" switches in the control histories shown in figure 2 were the result of high 
gains used in the damper logic. In addition, control surfaces were not rate-limited. 

This idealized damper system was used to avoid having to design a damper system for 
each configuration and to allow the same damper logic to be used for all configurations. 

Optimal recovery. - A six-degree-of -freedom optimization computer program was 
employed to investigate spin recovery. The objective was to obtain optimal control 
motions for rapid recovery from a spin and thereby to improve on the constant control 
results of figure 2. The computer program had a capability of iteratively determining 
optimal control histories so that some function of the state variables was minimized. 
Moreover, this minimization was achieved while constraints on specified state and con- 
trol variables were met. 

In the current application, the same spin conditions as those in figure 2 were used 
to initiate the recovery. Control variables were taken to be 6 e , 5 a , and 6 r . Con- 
straints were placed on these so that they could not exceed the limits given in table I. 

In addition, pitch-down elevator effectiveness and rudder effectiveness were set at zero 
for a > 50°. Terminal constraints were placed on a, \p, p, q, and r and were 
as follows: a = 10°, ^ = 0, p = q = r = 0. Thus, the aircraft was required to meet 
the same terminal conditions as in the recovery data of figure 2. Time was chosen as 
the performance index since, by minimizing time, recovery turns and altitude loss are 
also minimized. 

Optimization programs of this type require a nominal control history to initiate the 
iterative process. This nominal history is used to compute an initial time history after 
which the program adjusts the control variables on succeeding iterations in an attempt to 
obtain the optimal control motions. The nominal control history was one for which all 
control deflections were maintained at zero throughout the initial computation. 
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One additional quantity, a cut-off or stopping variable, is required by the optimi- 
zation program. The cut-off variable is normally taken to be one of the state variables 
such as angle of attack, yaw rate, or time. However, in the present application none of 
the state variables can be used to terminate an iterative cycle. Time cannot be used 
since it is the minimization function. Likewise, other variables such as a and r can- 
not be used since there is no assurance that they will reach the desired value (a = 10°, 
r = 0) during an iterative cycle. Hence, a "dummy" cut-off variable was adopted to 
terminate each cycle. This dummy variable corresponds to time on the nominal compu- 
tation, but on succeeding iterations the program adjusts the cut-off variable to achieve 
optimal performance. A cut-off value of 16 sec was used for the nominal calculation. 

After multiple iterative cycles, the final recovery results shown in figure 3 were 
obtained. As is shown, an impressive improvement in recovery time, altitude loss, and 
turns is achieved in that the constant control results were all reduced by greater than 
50 percent. 

Consider the optimal control histories shown in figure 3. The normal recovery 
mode of the aileron with the spin is used throughout the recovery. The rudder is used 
against the spin as soon as it becomes effective (a < 50°). The elevator history of 
figure 3 provides the most interesting study. Superimposed on the elevator history is 
the pitch rate 9 of the aircraft. During the initial part of the recovery the elevator is 
switched in phase with pitch rate. (It should be emphasized that the switching is based 
on 9 and not q.) That is, for negative 6, pitch-down elevator (normally 6 e ,butin 
this case, 6 e = 0°) was used, while for positive 9, pitch-up elevator was used. This 
procedure quickly excited an unstable oscillation in angle of attack which effects a rapid 
recovery as compared with the constant control technique of figure 2. (It is of interest 
to note that using the elevator to "rock" out of a spin was proposed in 1931 in ref. 9.) 

The fact that maximum rudder was not used when effective ( a < 50°) results from the 
slow terminal convergence characteristics of the optimization procedure. Had the iter- 
ative procedure been continued to obtain a closer minimum, maximum rudder deflection 
would probably have been employed, but little change in overall recovery performance 
would have resulted. 

In addition to the results of figure 3, optimal recoveries were computed for other 
aircraft configurations with various initial spin conditions. The control history results 
were in general agreement with those of figure 3; that is, the 9 and 6 g histories were 
similar to those of figure 3. 

Based on these optimal results, a suboptimal recovery logic was devised as follows: 

For a > a L : Rudder against spin 
Aileron with spin 
Elevator: Pitch up for 9 > 0 

Pitch down for 6 £ 0 
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For a ^ n L , the stability augmentation logic of table HI was used. Typical results 
obtained using this suboptimal recovery logic, referred to as pitch excitation, are given 
below. 

Pitch-excitation recovery. - A recovery using the suboptimal, pitch-excitation logic 
is shown in figure 4. Initial conditions for this case are the same spin conditions as were 
used on the constant control and optimal recoveries of figures 2 and 3. A comparison of 
figures 3 and 4 shows that the pitch- excitation recovery results compare favorably with 
those obtained optimally with only slight increases occurring in time, turns, and altitude 
loss. Thus, for the recoveries of figures 2 and 4, the pitch-excitation mode results in 
about a 50-percent reduction in turns, altitude loss, and time over that achieved using 
a constant control procedure. 


Configuration B 

Configuration B provides an interesting application of the pitch-excitation procedure. 
As was concluded in reference 4 and verified in the current study, constant control 
recovery from the spin conditions of table n is impossible without exceeding the control 
surface deflection limits of the aircraft. 

Constant control recovery. - An attempted constant control recovery history for 
configuration B is given in figure 5. Initial conditions for this case are those given in 
table n. Although no recovery was achieved, the control settings used yielded the 
greatest possible reduction in yaw rate and angle of attack. As shown in figure 5, the 
attempted recovery caused some reduction in yaw rate and forced the aircraft into a 
somewhat steeper spin. Extending the history of figure 5 to about 50 sec caused the air- 
craft to damp into another stable spin condition at an angle of attack and yaw rate of 
about 71° and 1.35 rad/sec, respectively. 

Pitch-excitation recovery .- A recovery from the same equilibrium spin conditions 
as in figure 5 was attempted by use of the pitch-excitation technique. The results appear 
in figure 6. As is shown, recovery was achieved after 24 sec and 5 turns. Thus, although 
no constant control recovery was possible for configuration B, pitch excitation effectively 
produced an unstable oscillation in angle of attack and returned the aircraft to a nonspin- 
ning condition. 

The recovery histories previously shown for configurations A and B are typical of 
those for other configurations and spin conditions investigated. Therefore, no further 
detailed time histories are included in the report. Results for other spin conditions and 
aircraft are, however, briefly discussed. 
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Summary of Spin Recovery Results 


Previously given recovery results for configurations A and B (figs. 2 to 6) are 
summarized in table IV. Additional recovery results are also provided for configura- 
tions A and C. 

A comparison is given in table IV between the constant control and pitch-excitation 
techniques for recovery from the run 2 spin condition of table II for configuration A. 
Although improvement in recovery performance from this flat spin condition using pitch 
excitation is not as great as for the intermediate spin, a significant reduction in turns, 
altitude loss, and time is seen to exist. 

An analysis of table II reveals that the flat spin characteristics for configuration C 
(run 4) are not as severe as for the other configurations; that is, the spin rate and angle 
of attack for the flat spin are consistently lower than those obtained with other models. 
These less severe spin characteristics for aircraft C are reflected in the recovery 
results given in table IV. Recovery is quickly achieved for both techniques, the pitch- 
excitation mode yielding some improvement over the constant control approach. 

Application of Spin Recovery Results 

The constant control, optimal, and pitch -excitation recovery histories presented 
should not be considered as representing typical or recommended piloting techniques. 
Although the constant control logic resembles recommended pilot procedures for initi- 
ating recovery from a spin during the early constant control portion, the subsequent 
switching in the damper logic is a highly idealized representation of piloting technique. 
Moreover, the pitch -excitation mode may not be a feasible pilot input even in the initial 
phase, since the large coupled oscillations would tend to disorient the pilot and make it 
difficult for him to respond accurately to pitch rate. Therefore, the suboptimal control 
logic should be viewed as indicating a possible approach to the design of an automatic aid 
to the pilot in recovering from a spin. 

In a practical system, it seems reasonable to assume that the pilot' s role would be 
to perform the final pullout after yaw rate had decreased sufficiently. However, precise 
definition of the pilot' s role and the form of a practical automatic-control system based 
on the pitch-excitation logic would require detailed piloted simulation studies and flight 
tests. Remotely piloted vehicles would seem to be ideally suited for initial studies of 
this type. 


CONCLUDING REMARKS 


Results have been presented from an analytical study of procedures for effecting 
recovery from steady spin conditions for three assumed aircraft configurations. The 
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performance of three recovery procedures which utilize conventional aerodynamic con- 
trols were analyzed with respect to altitude loss, turns, and time required to achieve 
recovery. Recovery procedures included a constant mode, an optimal mode, and a sub- 
optimal pitch-excitation logic patterned after optimal recovery results. Recovery per- 
formance using the suboptimal technique was shown to approach that obtained optimally 
and to yield a significant improvement over that achieved with the constant control 
approach. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 27, 1974. 
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APPENDIX 


EQUATIONS OF MOTIONS 


The six-degree-of-freedom nonlinear equations of motion used in the analysis are 
listed below. The aerodynamic coefficients are written in a form consistent with the data 
available for the most extensively tested model. 


Normal force: 


v . /ton as C ZM) + C Z 6 M 6e 

O' = q + (tan a tan /3)/3 - ( a W + -2— 

\ V ) mV cos a cos ft 


gfcosecos 4> 


V\cos a cos 0 


tan 0 
cos ot 


Lateral force: 

^ = -(^)v + —^—^IC Y (a,0) + Cy^M) 6 e + CY 6a (a,/3) 6 a 


qS 


+ C Y 6r M 6 r + ^[CYp(«) P + c Y r («) r] + f ^ 


+ p sin a - r cos a 


Longitudinal force: 


« 

V = 


(V tan a) a + (V tan 0)fi + 


q S c x( a >® + c Xg e (° ! ^) 6 e 
in cos a cos /3 


g sin 9 
cos a cos |8 


qV tan a + rV 

cos a 


Pitching moment: 

4 = (^I7^) Pr + ^ r2 " p2 ) + + C m 6e («,0) 6 e + ^ C mq (o!) qj 

Rolling moment: 



+ C l6r (a,(3) 8 r + ^[C/ p (a) p + C lr (a) rjj 
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APPENDIX - Concluded 


Yawing moment: 

» - (f)e - ($)V - * f (c»<M> ♦ Cb* S e + C^<*» a a 

+ C n 6r ( Q! ,/3) 6 r + ^[ c n p (a) p + Cn r (a) r j 
In addition, the following formulas were used: 

6 = q cos 0 - r sin 0 

0 = p + tan 0(q sin 0 + r cos 0) 

• q sin 0 + r cos 0 
cos 6 

Turns = \ 0 dt 

u = V cos a cos /3 
v = V sin /3 
w = V sin a cos /3 

h = u sin 9 - v sin 0 cos 9 - w cos 0 cos 0 
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TABLE L- WEIGHT, INERTIA, AND GEOMETRIC CHARACTERISTICS 



76.247 4.465 - 13.963 3.513 86.992 .271 9 144 

73.210 - 2.383 - 16.810 -.365 78.550 .188 9 144 




















TABLE III.- STABILITY AUGMENTATION LOGIC 


3 

VII 

a 

6 r = lOOOr 


6 a = lOOOp 

r > r L : 

5^ = lOOOq 

r = r L : 

6 e = 5 fa - a-p') + 100a 


TABLE IV.- SUMMARY OF SPIN RECOVERY RESULTS 


Configuration 

Run 

Recovery technique 

Turns 

Altitude loss, 
m 

Time, 

sec 

A 

1 

Constant controls 

8.3 

3250 

28 



Pitch excitation 

4.0 

1750 

14 



Optimal 

3.8 

1550 

13 

A 

2 

Constant controls 

12.9 

4000 

36 



Pitch excitation 

7.3 

2220 

19 

B 

3 

Constant controls 

— 

No recovery 

— 



Pitch excitation 

5.0 

2050 

24 

C 

4 

Constant controls 

3.3 

1800 

22 



Pitch excitation 

2.6 

1400 

17 


17 





18 
















13 , 

deg 


deg 


deg 


P; 

rad/ sec 


q> 

rad/sec 


r, 

rad/sec 


30 

0 

-30 

60 

0 

-60 

0 

-40 

-80 

3 

0 

-3 

1 

0 

-1 

0 

-1 

-2 
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